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I. INJECRUCTION 


A variety of acoustic models exist to predict transzis- 
sion less. Each of these mcdels contain inherent strengths 
and weaknesses. All have shown poor results in a shallow 
water envircnoment due to difficulties at the ocean Erctton 
where there is an interface lLetween media of differert sound 
speeds ard densities. 

Since its introduction (Hardin and Tappert,1973), the 
parakclic wave eguatica has been a widely accepted szeans cf 
soluticn for accustic propagation. The earliest frecrams 
used a sgplit-step Fcurier transform alyorithni to sclve the 
Farakclic eguation (FE). Severai other solution technijues 
have reer develored ;rimarily to overcome difficulties that 
cccur sken the Fourier transform encounters an interface 
ketween different media (Lee ard Botseas, 1982 and ¥cDaniel 
and lee, 1S&2). 

An alternative sclution technique that uses an iuplicit 
finite-difference (IEFI) algorithm was developed by lee and 
EFapadakis (1979). Ibis methcd incorporates apfrogriate 
interface conditions and allows solutions in shallow water. 
Starting with the IFI algoriths, Jaeger (1983) develored a 
computer acdel to yredict transmission loss and accustic 
Fresstre tased cn user specified bottom tofography ard a 
single scund speed profile. This computer model uses the 
mathexatical treatmert of the horizontal and slopirg irter- 
faces developed by #&cLaniel and Lee (1982) aco lee and 
McDaniel (1583). It also utilizes several design features 
and wmetkcdse inccrporated in an earlier conputer fprcgran 
develczjed Ly Lee and Eoteas (1682), and a PE computer sodel 
develcgped ky Breck (1578). 
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Ike IFC program preserves continuity of pressure and 
contiruity cf the acrmal ccmpcnent of particle velccity at 
an interface between nedia having different scund speeds and 
densities. This feature makes the prcgram uaconditicrally 
stable and better aktle to handle tne bottom Ecurdary 
conditicr. 

Since its deéevelcjment, the IFD grogram has rect ceen 
rigorcusly tested. This thesis analyzes the jregrag's 
fecfcruarce in an idealized shallow water environment. The 
envircrgeert includes a simple siofing saad bottom breneathé ar 
isosreed water field. The aralysis begins by comparing the 
IFD‘'s Eredicticns with predictions frog Jerscn and 
Kuperugan's (1980) FE model and Coppens, Humphries and 
Sander's (1684) imace model. The analysis also includes a 
compariscn of the acdelts estimated transmissicn icss 
contcurs with expectations based On sipie Ehysical 
reascning. Finally, the thesis describes an attempt tc scdei 
a shallcw sloping bottom in an experimental tank. The tank 
contains a sand bottcm sculptured with a ten degree sloje. 
laboratcry measurements ci tae pressure field are taken at a 
freguency cf 100 kHz for ccmgarison to the predicted rfres- 
sure field cenerated by the IFD. These comyariscns cf the 
IFD predictions witk cther acdel estinates, theory, and 
laboratory measuremerts yive an indication of the IFLI's 
perfcrmance in a Shallcw water environment. 
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II. JHE IEE CCHEOTER SODEL 


a. =6EACKCECOND 


Ac igglicit finite-difference solution technigue tc the 
Farakclic eyuation has been studied and refined ty many 
authors. The history of this development is explained in 
detail by Jaeger (1585), but merits review in order tc yain 
a perspective on an aralysis of the IFD computer model. iuae 
Farakclic eguation is an ayproximation to the eliirftical 
wave é€guaticn. The first means of solving the PE useda 
split-step fast Fourier transform method as develcred by 
Tappert and Hardin (1973). This method reguires periodic 
Eoundary conditicns ix depth kecause of the finite Fourier 
transfcrm and handles this constraint Ly introducing ar 
artificial horizontal fressure release bottom belcw tae 
actual fhysical bottcn. This method of implementing an 
artificial bottom was inccrgcrated into the earliest PE 
models developed by Jenson and Krol (1975) and Brock (1576). 

Errers in this sflit-step Fcurier transform methcd were 
found tc be progorticnal tc the horizontal range ster ana 
the second derivative of the index of refraction. The séeccnd 
derivative of the izrdex of refraction tends to te large 
across tké€ ocean bottcm interface. Another problem with the 
split-step method is that it does not consider density 
differences between two different media at an interface, 
which influences the reflection coetficient. For these 
reascnse the split~ster fourier transform methcd preved to be 
Foorly suited for a skallow water environment. 


Ibe 1IFL soluticr metbod was introduced by Lee and 
Fapadakis (1979) as ar alternative to the split-ster methcd. 
The IFD method emplcys a seccnd order central difference 


i forzula to solve the FE in the fcrm of a tridiagonal matrix. 

7 Althctgh the first version of the IFD did handle discontiru- 

ities in the scund speed profile, it did not consider 

density discontinuities. In 19&2 McDaniel and Lee intreduced 

a metked tc hardie a horizcntai interface of different . 
densities. In 1983, they extended their treatzent to 

include a sloping interface. It is this version of tne IFD 

: that is used in Jaeger's computer proyran. 
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En. TEE CCBEUTER MODEL 


The JFL computer yrogram consists of a main prograz end 
= twenty sukrcutines. the pregaz utilizes a modular ccnstruc- 
e tion sc that each of the various subroutines are called from 
the main program to ccaupliete a specific calculation cr furc- 
tion when reyuired. The IFD is run interactively fro a 
oy user generated ingut file that contains values for 
= frequency, one sourd speed profile, a bottom jfrecfile, 
sourceyreceiver depths, attenuation coetficients for Ecth 
the water and the bottem, and several other input farameters 
e that tell tke the prcecram where to obtain a solution within 
S the field. The procram initiates the calculations assuming 
an initial Gaussian jyressure field and an artificial ;fres- 
sure release surface at a user specified depta. 

Atteruation in toth the water and tae sedisrert is 
handled as complex indices cf refraction. An artificial 
attentation layer is estaklished beneath tbe sedixert to 
introcuce attenuatior above the artificial pressure release 
surface. The actual gagnitude or this enhanced attenuaticna 
is calculated using ar equation derived by Brock (1978) for 
use in his FE computer nodel. 

Ike IFD program steps along the specified bottca ;rcfile 
and sakes calculaticns down tnrouyh the water/sediment 
column at €ach user syecified horizonal ranje. The prcgraa 
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requires that the Ecttcm intersect exactly at a vertical 
grid reirt. As a result, for a Sloping bottom the grejraz 


- automatically calculates the range step to fulfill this is 
reguirement. This ccmputer generated range step is always ae 
= ; less than cr equal tc the user provided range ster. As the ns 
ee Slope cf the bottom increases, the range ster must decrease ae 


and sore calculations and computer time are reguireé to 
solve tke entire pressure field. In scmne Situations an 
actual mcdification cf the user inputted bottom is reguired. 
Tnis cccurs with a very gently siopiny bcttom when the ee 
reguired rarge step exceeds tbe user specified range step. ae 
Here, the programs automaticall; modeis the bottcm as a ‘ 
series cf level and sloping sections in order to ensure the 
user cenerated range step is net exceeded. This mcdification 
cf tke Ecttcm is always less than or egual to one-half the 
vertical grid spacing and the podel issues a warnirg tc the 
user cf the modificaticn. 
Ecth printed and graphical output are provided ty the 
p= IFD. The Frinted output prevides transmissicn loss and the 
real and ifaginary ccaponents cf the pressure field at each 
depth fcr a specified bcrizcaotal range. The graphical cutjut 
is a jlct cf transzission loss versus range at the user 
specified receiver depth. 


Cw. &&CLEL EROBLESS/MCLIFICATICNS 


As this study of the IFD pecgram progressed, it became on 
necessary tc modify certain aspects of the model. Mcst of asia 
: these scdifications sere necessary to alter the fregrana - 
n cutput irtc a mcre desirable fcrn, but a few were insple- 
gented tc correct picgramming deficiencies. Although this 
F sectican cf the thesis discusses the earliest mode] rurs, Toe 
iS these results are rct presented in detail, but orly 
cao discussed in general terms kecause they were obtained tefore 
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the ccrgfuter model tas fully accdified. All modificaticns 
sere ace coaly after careful analysis of multiple aodel 
runs. It is impertart to realize that all results generated 
for ccrparison to ctker models and Laboratory measuresents 
came frcz a fully mno¢cified version of the IFD. 

Since the ultiszate otjective of this study was to 
compare 1FD model results with laboratory measurements the 
Frogras sas first cur with irfput parameters which exactly 
modeled ccnditicrs in tne tank. The experimental set uf, 
explained later in great detail, consisted of a tazk that is 
appreximately twe meters in length, one meter in depth, with 
aten ¢ecree slcping sand bottcm and a paaxinum water derth 
cf 35 centizeters(cm). Based on a test case run ty Jaecser 
(1983) cf a simple sloping bcttom, and by results shcun by 
Jenscr ard Kupergzan (1980) for propagation in a wedge-shared 
ccean, it was expected that there would be certain reéecccniz- 
able patterns in the predicted Fropagation fatterns. 
Specifically, since tke speed cf sound in the bottom exceeds 
that in the water (fast bottcs), the siafle sloping Ecttoa 
supperts trapped nacraal scde propagation (Coprens and 
Sanders, 1981). As the acoustic energy travels upesicpe 
toward tke apex, successively lower modes are cut cfi and 
the enercy contained in these sodes is transmitted intc the 
kottcr. The range from the apex at which enersy cf the 
lowest scde is transsitted intc the bottom is referred to as 
the duar distance anc is a function of wavelength, wedge 
angle, and the ratic cf sound speed in the water tc the 
sound sreed in the sediment. An empirical equaticn ttkat 
defines this dumg distance was derived by Coprens, Sanders, 
Ioanncu, ard Kawamura (1978); and was used to identify the 
expected ranges of these durp distances for the sciven 
scenario described aicve. 

Ike initial unmodified IFD run used the parameters taken 
from the tank and showed ro recognizable patterns ir the 
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acoustic field. There was no observable decrease in trars- 
tissicn lceS at the various dump distances as exyected. 
Rather, results indicated widely fluctuating patterre in the 
acoustic tield that appeared inconsistent with both rrevicus 
studies and sizuple physical reasoniny. Upon closer analysis, 
it was disccvered that althcugh the prograz was designed to 
te independent of scale there are several logic statezents 


tet atet a a 
e 


that are not implemerted if the user provided range steyz is 
, less thar cne meter. Because the logic statements aren't 
satisfied, the NEWSEG and MEWMAT subroutines (Jaeger, 19&3) 
are nct called cerrectly. ‘The NEWSEG subroutine initializes 
the Ecttcm slope anc the NEWMAT subroutine computes matrix 
€lements for the fprecraa. Obviously, errors in these two 
Frogrags functions sericusly distort results. Because of this 
systematic error in tke progragz it became necessary tc scale 
up all tank parameters. Ail distances were scaled uf ry a 
factcr of 1000, and frequency was scaled down by a gactcr of 
a 1000. Careful analysis reveals that all input parameters are 
- a functicn cf either distance cr frequency, so this scaling 
Froduces results that sodel thcse expected in the tank. 

Ike second modification cf the IFD was reguireé to 
Frovide a three dimersional grarhics display. As discussed 
earlier, the IFC prevides a transmission loss flot versus 
range ata single depth. Hcwever, to study the srodei 
Fredicticnse in greater detail it was felt that a twe dizen- 
sional aralysis cf tke model estimates would be acre sean- 
dngful. As a result a transmission loss contouring frecgranm 
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was developed. The yrogram (Appendix B) displays transmis- 
sion lcss contours for range versus depth. Use cf the 
contcur sleet reguires that transmission loss values gener- 
ated ky the IFD be sent to a data disk used by the contcur 
routire. Ic facilitate this transfer a dummy variakle (LifD) 
was estallished in the PRINI2 subroutine to store the trans- 
sissicn loss values and then these values are writter tc the 
Gata disk at the end cf the gain frograag. 
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Anctker rodification of the IFD output was reguired to 
Chance the real and izaginary: components of the pressure to 
asincle fressure auzlitude sagnitude. It was felt that 
Gealing with the ypressure magnitude was easier ard scre 
: geanincful than witk the real and imaginary components of 
: the pressure field. This ccnversion was done in the PRINT2 
subroutire and established a new variable (PRMAG) tc repre- 
sent tke pressure macritude. 


A final modification in the computer prcgram was made 
. due tc a suspected errer in the computation of the attenua- nee 
tion in the artificial layer. Fhysicai reasoning dictates re 
that preper iaplemertation of the artificial attentuation 
would result in a steady drcp off in acoustic pressure with 
deptk thrcushout the artificial layer, with pressure édrop- 
Fing tc zero at the fressure release surface. IFC accel 
results cn the other -and showed wide fluctuations in fres- 
sure sith depth in tke layer and then only a ainizal fall 
cff at tke pressure release surface. The eguaticer in the 
IFD that actually cozgutes the saynitude of the attenuation 
in tke artificial layer was taken directly from Ercck's 
i (1978) FE acdel (Jaecer, 1583). However, Brock's esguation 
was derived with feet as the unit of measurenent while 
daeger's ucdel is deéezived witb meters as the unit cf seas- 
urement. With this in nind, Jaeger's eguation stcula ke 
apprcxzimately a factcr of three larger than Brock's eguatior ‘ 
to ccrrect for the difference in units. To correct for this Pas 
error tke ecuaticn tc calculate attenuation (ATT(I)) in the ae 
NEWMAT sukrcutine was increased by a factcr of three. wien 
this ccrrection was isplemented, the large fluctuaticrs ina 
Fresetre with derth were elipinated. The expected dror cff ie 
in pressure with depth and the fall off of pressure tc zero poe 
at the pressure release surface were noted. 
A listirg of the revised IFD computer fprograa sith all 
sodificaticrs can be seen in Appendix A. 
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IL. MCLEL VERIFICATICS 
1. Cegpariscn with Jaeger Model Run 


In light of the nodifications to the IFD ccnryputer 


Frograg just discussed, it was necessary to ensure that the 


the chances thesselves did not introduce errors irte the 
wodel. ‘£0 asa first step the modified IFD rogram was run 
for cne of the test cases used by Jaeger in his originai 
work. Tkis case analyzes propagation in ar environment that 
goves frcm deep to skallicw water. This particular envircno- 
sent is depicted in Figure z.1 and was chosen because it was 
very Sinsilar to the simple slosing bottom in the tank exyper- 
iment. A sclution is obtained for a bottom with an ufsicpe 
cf 8.5 degrees. An isospeed water field is used witi scund ; 
speed set at 1500 ays. Scurcée frequency is 25 Hz. ihe e, 
scenario has a maximus depth of 350 meters and a range cf 40 
kilcweters. Both the source and the receiver are set at 25 
eters. hee 
The results using the rodified IFD program were the Pee 
same as Jaeger's for this test case above the artificial a 
atteruation layer. As explained earlier, the prcgrags was eee 
gsodified tc reflect higher attenuation in the artificial cee 
dayer, £0 as to pregerly reflect the effect of attenuation 
in this region. The sodified program results did sect show ae 


the large fluctuaticrse in pressure with depth in the artiti- ere 
cial layer but ratker the gradual decline in fressure ey 
towards a value of zero at the pressure release surface. Eut 


above tke artificial layer the modified program results were 
identical with those achieved ty Jaeger with the criginal 
Frograa. Apparently, the sincr caanges in the pregraa 
designed tc improve cn the form of the program output dces 
not binder the acdelts ability to achieve a soluticr in the bee 
upper sediment and water column. ve 
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FREQUENCY - 2.0 HZ 
SOURCE DEPTH - 25 M. 


SEDIMENT 


ARTIFICIAL LAYER 


20.0 
RANGE(km) 


Figure 2.1 Jaecer's Deep-to-Shallow Water Case. 


ensen and Kuperman Model Run 


Ihe seccnod attempt at verifing the IFL scadel 
involved a compariscr of mcdel results with those achieved 
with a FE wcdel desicred by Jensen and Kuperman (1980). The 
Jenséxr and Kuperman Mcdel (JKM) uses a split-stefr solution 
techrigue. Compariscr with this particular model was chosen 
kecause it is one of the few that obtains a soluticr ir two 
Gimensicns. Most acotstic acdels obtain a solution at only a 


Single depth. In addition, Jensen and Kuperman sade their 
godel rurs in a simple sloping ccean bottom envircrment very 
Similar to the scenario of interest sodeied in the tark. 
This e€rvircoment is depicted in Figure 2.2 and features a 
gently sloging bettcs of 2.2 degrees. The water coluzn has a 


eyo 


unifcrs speed of 150C sys. The source is placed just Lelow 
the aidgcint in the channel at 112 meters and is driven at a 


rf frequency cf 25 Hz. The maxinum depth in this scerario is ; 
z00 seters and has a gaximug range of 12.5 kilometers. wee 
. Ihe Jensen ard Kuperman results were expressed as See 
FREQUENCY - 2.5 HZ. mo 
SOURCE DEPTH -112 M. on 
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: oe 
: Pigure 2.2 Jenser and Kuperman Slofing Bottoa Case. ae 
: a 
- transmission loss ccrtours fcr range versus depth. MTheir 
5 study ccncentrated cn transrission loss patterrs in the Sao 
: sedizent, rut results were cbtained both in the sedizent and ae 
in the water. These results are compared to the IFD results one 
. at ranges cf 2.5, 5.C€, 7.5, 10.0, and 12.5 kilometers. ‘These — 
: ranges were chosen for analysis because the JKM results oe 
. showed the greatest variaticn in transmission loss with ee 
: ; eed 
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depth and thus, sakes for a more meaningful compariscn with 
the IFI results. 
The IFD and JKM results can be seen in Figure z.3 
| through Figure 2.7. In all the figures, the IFD estizates oe 
are skcwr as a solid curve thile the Jensen and Kupergan : pee 
results are depicted as circular points. The first analysis ibe 
is at a rance of 2.5 kilometers. From Figure 2.2, it can be oe 
seen that tte depth at this ranye is 200 meters and is ira ais 
flat kottom reygion. From the results shown in Figure 2.3, 
it is ckvicus that [cth mcdels obtained almost identical 
results frcem the strface dcun to a depth of abcut 300 - 
meters. There are differences between the two sets of 
Fredicticns from the ccean kottom to a depth of 100 neters 
telow this ysoint. EFelow 300 seters the Jensen and Kurperagan 
: results shos a very slight iacrease in transmissicn Ilcss 


L (IL) witk depth. Tke IFD alse shows an overall increase in ere 
: transsissicr loss with depth tut with several fluctuaticns pees 
in trarsgission loss and a gsarked peak at about 300 néters. oon 
; So in general, the results frcm the two models have the sane _ = 
i generai terdencies altkough the IfD appears to shcw greater 
: detail in results near the water/sediment boundary. : 
Figure 2.4 stcws results at 5.0 ka in range. Here, 
7 the water depth is still 300 meters and sarks the very 
; teginrirg cf the slozing bottos section. For this range the Sos 
; JKM predictions are crly available to a depth of 300 seters. mae 
: The results are nearly identical with those cbhtained ry the 
; IFD. Bcth models show a relative minigua in transaission ee 
4 loss at a depth of 1C€0 meters and then a gentle increase in re 
; IL with depth. 
: The aodel results at a range of 7.5 km are seen in Bese 
; Figuce 2.5. At this range the bottom depth is akcut 150 as 
P meters ané the bottca is slcfpiny. Ine models show the aa! 
; greatest difference at this range. Both gpodels freduce oe 
: nearly identical results in the water coluan, but bkeneath verte 
: st 
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Pigure 2.3 IFC and JKM Ccaparison at a Range of 2.5 Ka. 
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Pigure 2.4 IPD and JKS Conparison at a Range of 5.0 Ka. 
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Figure 2.5 IFC and JKM Cemparison at a Range of 7.5 Km. 
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Figure 2.6 IFC and JK Comparison at a Range of 1C.0 Ku. > 
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150 setters in depth the mcdel estimates ceyin tec stew 


Cifferences. The Jensen and Kuperman results show a gradual “ 
increase ir TL with depth below the ocean kottoa. The 1FD Sa 
results cn the other tand, shcw a gradual decrease in trans- 
kissicn Ices frem 15C neters to 300 meters, with a relative 
Qinigusa at 300 seters and then an increase in IL reneath 


in tke sedineat, hcwever, beth models produce sisrilar 


this de;th. Althouck there are some differences in scdel e 
estizates near the kcundary in the flat pottom secticns 7 
previcusly discussed, the differences appear to ke greater Se 
at the tourdary in this slcfing bottom case. Ir _ thcse ee 
regicns away from the Ecundary in either the water cclugzn or ae 


results. The JKS's difficulty in obtaining an accurate solu- 
tion at tke sedizent Ecundary is not totally unexpected. ithe 
model cses a variaticn of the straight split-stef solution 
technigue similar tc Brock*s computer model (Jenser and 
Kupersan, 1980), tkrat has characteristically been unsuc- 
cessful in cbhtaining a reliable sclution near a boundary. 
The results at a range of 10.0 km are seen in Figure 
226. At this range the water depth is apfroxisately 100 
weters and the bottog is again in a sloping region. At this 
range Jensen and Ktyerman results are only available to 
about 15C seters in depth. For the data availatie, the 
models prodtce nearly identical results. Both models predict 
@ minizug at about £C meters in depth and then an alacst 
linear increase in Tl with depth. The IFD results also 
reflect 2 transmissicr loss sazigum at about 200 meters, and 
then a slight decrease of IL beneath this depth. Data at 
these depths are not availakle from the Jensen and Kuyeérgzan 
run. In dJight of tke results cbtained at 7.5 ka for a 
sloping Ecttom case cne might expect that the resylts fer 
the two scdels weuld show differences near the ocean t-ottca. 
However, since there is crly one Jensen and Kupferman 
Fredicticn available near the fottom for this range, it is 
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difficult tc aake any definite conclusicns regardizg ditfer- 


ences ix acdel perfcraance at the boundary. 

Results at tle apex (range equal to 12.5 kz) can be 
seen in Figure 2.7. Since all estimates are nade in the 
sedigent at this rarge, there is no boundary tec ccrtend 
with. At this range the two models snow the test agreenert. 
Eoth models show a gradual decrease in transmission lcss to 
a derfth cf about 300 meters and then a gradual increase cor 
IL with depth. 

In general, there afrears to be good agreement 
ketween tke IFD and the JKM sodel results in regicrs not 
influenced ty a bounéary. In koth the water and the deep 
sediment the twce models preduce Similar results. Although 
this is ret conclusive evidence, these similarities sucgest 
that the IFD can successfully model accustic propagaticn in 
these recicrs. Near tke water/sediment boundary however, the 
Jenser and Kupergan ard IFD predictions show marked differ- 
ences. These differences afrpear to intensify as the iEctton 
keccueées mcre sofhisticated. In general, the Jenser and 
Kupeérgan results dc rect seem to show the detail the IFD 
results do. Considering the different solution technigues 
emplcyed Ey the two acdels, these differences in results are 
expected. 


3. Cegrariscn with Coprens,Humphries and Sarders Rodel 


im 
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the third attempt at verifying the IFD Lperfcruance 
was dcre ky comparing results with an image theory srodei 
derived ty Coppens, Humphries, and Sanders (1984). This 
Coppens, Humphries, and Sanders Model (CHSM) uses a saddie 
Foint agrreximation to an image model te solve fcr the 
acoustic field. Bcth programs were run for the scaled 
scenaric mcdeled in the tank experiment. This scenario is 
depicted in Figure 2.8 and features a ten degree slcping 


2é 


ktottcz. The maxignug water defth for the run is 350 séters 
and the saximum range is twe kilometers. The source gener- 
ates a 100 bz signal and is set at a depth of 175 meters. 


FREQUENCY -100 HZ. 
SOURCE DEPTH -175 M. 


SEDIMENT a 


0.5 1.9 
RANG=(km) , 
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Figure 2.8 10 Legree Sigsple Sloping Bottom Case. 


The image model provides soluticns only at the apex. 
The ccupariscn of fredicticns at the apex(range egual tc 2 
kilozeters), can be seen in Figure 2.9. In the Figure, the 
IFD values are slotted as the solid curve while the image 
sodel values are shown as circular points. The predicticns 


are stcwr as values cf norgzalized pressure amplitude fora 
given dezth. The difterent models produce pressure agplitude 
values in different writs, and thus, had to be normalized to 
gake a ccmgarison of values possible. This normalization was 
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Gone for €ach model ty dividiny tie presssure arplittce by 
the maxizugz amplitude predicted by the model. 

fhe results suggest that even thoujh the scdéeis 
predict sisilar larce scale trends in fressur2 auplitude 
with depth, there are several differences in detail. The 
Flot ¢hcus that toth the IFL and the igzage mcdel predict aa 
almost linear increase in auafplitude with depth Coun tca 
specific saximum and then a Slow decrease in ausglitude 
reneatl this maximum. The IFD however, shows tae gaximum at 
about 15.5 meters in depth while the imayeé model saxigun is 
deeper at 1S meters. Beneath this naximunm the IFD shcws a 
wucn skarper decrease in amplitude than the image rcdel. lhe 
trend in the image mcdel data appears smooth, wnile the IFD 
curve reflects several szall scale fluctuations. 

It appears tkat both models predict sinailar large 
scale trends in presstre amplitude for this scenario. It is 
difficult, if not imgcssible, to account for the differences 
in detail. As a minigun, hcwever, this comparison indicates 
that the two sets cf predicticns are consistent witk cne 
another and can be considered reasonable in this shailow 
water envirconent. 


4. Ccogsariscn with Physical Keasoning 


Frcun model ccrparisons it appears that the iIFE at 
least makes a reasonalle estimation of acoustic fields ira 
shallew water envirorzent. IFD results are also analyzed in 
comparisecn sith basic physical reasoning and theory as a 
fourth attearpt at gwodel verification. This analysis ¢exazines 
IFD transzission loss contcurs for the simple ten degree 
sloping ccean scenaric seen in Figure 2.8. The scenario is 
two kilczeters in rance with @ maximum depth of 350 meters. 
The tke source generates a Signal at 100 Hz and is placed at 
175 metere in degth. 
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Figure 2.9 IFC and CHSM Comparison at the Apex. 
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Irarsmissicn loss ccntour plots are diszlayeeé in 
Figure 2.10 througk Fiyure 2.13. All TL contours are 
expressed in decibels(db) and are shcewn as a functicn of 
depth wersus range. In ali figures, the ten degree slcping 
kottca is shown as a solid unlabeled line. The contcur 
Flots are displayed over fcur range subsections due to 
computer craphics ligitaticns and to emphasize different 
significant features in the field. 

the first figure shows transmission ioss ccrtotrs 
from the scurce to a range of 600 meters. Ine ccntcure is 
the first 250 meters appear very symmetric, increasinj; 
Cutward from the sotrce in a yattern that resenbles spher- 
ical sepreading. Since the IFIZ program assumes a Gaussian 
starting field this early pattern is expected. 

Freg a range of about 300 meters to 60) meters the 
field in the water appears tc be dominated by a surface 
reflecticn pattern. From the spotty appearance cf the TL 
contcurs in this regicn there is an indicaticn of ar irter- 
acticn cf surface reflecticn and bottom reflecticn cn the 
contcurs. It is sossitle to ccapare transmission loss uaxina 
with rodes in the surface interference pattern. Easéd on 
surface interference theory these nodes should cccur where 
(Kinsler, Frey, Coppers, and Sanders, 1982): 


SIN(khdyr) = 0 
ex in cther words; 
khdsr = nf 


where; 

pB=0),1,2 « « « 

r= Range 

= Wavenumber (2,9) 
Scurce Depth 

= Depth cf Node 

= Wavelength. 
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Ey maripulating the akove €,uation it is ypossible tc sclve 
for tre deyth at which the acdes should be observed. for 
this particular scenario at a ranye of 604 meters, rodes 
should cccur at integer multiples of 25 meters in éefth 
(25n). Figure 2.1C shows that these transaissicn icss 
maxima dc exist as expected, every 25 meters in depth at the 
stated range. 

Fased on simple physical reasoning one weuld expect 
refractive Lending aleng the ocean bottom due to density 
differences between tke water and the sediment. The certcur 
Flot reflects a charge in tasic pattern at the interface. 
There is a Lending of the T1 ccntours that suggest a refrac- 
tive influerce. 

Figure 2.11 shows transmission loss contours frog a 
Tange of 60C meters tc 1350 mé€ters. Again in this recicr the 
water aprears to be deminated ty surface reflection. Sclving 
for the nedes in this surface interference pattern at a 
range of 1350 meters, it is found that taese nodes stculd 
appear every 55 wetere in depth. From the figure it is again 
seen that Tl maximuas do occur every 55 meters in degth as 
anticipated. As in the first figure there is a change ix 
the kasic IL pattern at the ocean tLottom. The rerding 
appears score accentuated thar in the previous figure, tut 
still suggests the influence of refraction at the cceéan 
kottca. 

Figure z.12 displays transmission loss frem 1:50 
meters tc fast the apex at a range of 2100 meters. In this 
figure, tke dominance cf the surface reflecticn mechanism is 
Jess ckvicus and the Il patterns become more complicated. 1t 
is ik this region that the influence of trapped ncerzal acde 
Fropagaticn can te seen. As discussed earlier, as accustic 
energy travels up tte slope tcward the apex, norzal sodes 
are cut off and energy is transmitted into the [cttca. 
According tc adiabatic normal mode theory, mcdal separation 
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: is rarge deyendent (Graves, Nacsel, Uberall, and Zaur, 1£75 oa 

= and Ccprens and Sanders, 1580). The ranye from the apex at . ee 
which tke lcwest ucde is transmitted into the bottom, can be 
e calculated using the fcllowing equation (Coprens, Sarders, 
: Icanncu, and Kawamura, 1978): 


= Aya SING TANG 


wiere: 
Dump Distauce Of The Lowest Mcde 


5 = Wavelergth 
6 = Critical Angle 
@ = Wedge Angle. 


Accoréing tc adiabatic norgzal mode tneory (Kinsler, Frey, 
Coppers, and Sanders, 1982), in deep water (near the scurce) 
the lcwer normal modes are far above cutoff and the adia- we 
katic eigenfunctions consist cf am inteyer number cf nalf oo 


: sine waves with zerc fressure at both the top and fEctton 
: surfaces. At the cutcff of each mode, the pressure at the 
a tottcz sust be gaximized; resulting in an adiabatic eigen- es 
; functicn that contains 1/4, 3y4, and 5/4 wavelengths at the ae 
respective cutoff distances cf 1X, 3X, and 5X for the three or 
dowest wcdes. From fFigure 2.14 it can be seen that as the 
nermal mcdes travel ug the slcping bottom, succeéesive modes 
‘ are fcerced into the Icttom at distances where a particular 
mode reaches a deptk at which it can not longer frecfacate. 
= Also frcm the figure it is chvious that a source set at 

gid-déepth can net excite the second node. Since tkis 


Farticular ceometry is present in the tank scenaric it is oe 
expected that the energy asscciated with this seccne ocde 


should rct be seer. Based on this line of fEhnysicai 


Teascring, modes shculd be dumped into the sediment at the Tare 
re first dump distance, fifth dump distance, ninth dump na 
ey distarce and so cn. If these modes are dumped as descrited a 
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Figure 2.14 Normal Bode Prorogation in a Wedge Shared Ocean. >... 
then they should be seen in the contour flot as reducticns oe 
cf trarcgxission loss in the bottom at ranges of 1952 meters, ey 
1760 seters, 1568 meters, ard so on. oe 


Frea Pigure 2.12 it can be seen that there is an 
cbvicts decrease in transmissicn loss at approximately 1£50 
meters (dump distance #1) and 1760 meters (dung distarce 
#5). Ibere is also a less clearly defined reducticn in IL > 
along tre tottom at 1565 meters (dump distance #9). 

The final cortcur plot (Figure 2.13) shows transmis- 
sion Icss contours from 1550 meters to 2300 meters. This 


figure is ar extensicr of Figure 2.12, intended to erphasize Rs 

how clearly the beam at the first dump distance is defined Fr 

in the slots. Ine kéam at the r£ifth dump distance ig also ; oe 
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visitle rut is not nearly as well defined. Fren this figure, 


“. ene car als¢c see an indication of a very narrow keagz ir the 
sediment at about 18€C meters. This distance corresrcnds to 
the trird detmp distance (second normal mode). Based cr adia- 
katic scde theory tke seccrd mode is not expected tc te 
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excited. Hcewever, adiabatic mode theory is only an aprroxi- 
maticr cf rcermal mode behavior. This approximation cf nerasal i 
pode tketavior teccuzes less exact as the bottcs slcre we 
increases and the clcser the source is to the ajex. the MELE 


appearance cf a narrcw beam at the third dump distance indi- 
cates that there is a strceng fossibility that the seccnd 
: acde is present and tkat the adiabatic approximaticr is ret 
oe exact with a ten degree bottcz slope. 

5 fhe basic features of the IFD contour flicts are 
consistent with botk physical reasoning and theory. £Easic 
surface reflection and tottom refraction occur sséere 


so 
eee 


a expected ard behave as anticipated. In the far field, 
trapred norgzal mcde fropagaticr is observed and can te veri- 
fied with sigple dumg distance calculations. The locaticn of 
Teams dusped into tke bottom appear consistent with tasic 
rs : mode thecry. In short, the transmission loss contours indi- nai 
. cate that the IFD is making reasonable predictions cf the Yer 

acoustic field in a ¢kallow water environment. — 


Z 5. Verification Summary oes 


It is difficult to say how exact the IFD precicticns 
are fcr a shallcw water environment based on these single 
_ verificatior technigues. As a minimum it can be said tkat 
. the acdel] results are at least consistent with other godel en 
5 Fredicticns and exgectations tased on simple ghysical hee 

reascnring. Model estimates are virtually the same as the 

Jenser and Kupergan FE model in regions not influenced cry a 
os waterysedigent toundary. Clcse to the boundary the IFD hats 
2 results appear to shew greater detail and variaticn than an 
"| 
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; 
‘@ 
e this EE rcdel. The JFD results aiso appear consistent with 
. the general trends in ;ressure agplitude predictecd [Fy the 
7 Coppens, Hugphries, and Sanders imaye model. Agair differ- 
a ences were roted in tke small scale structure. Finally, the 
‘ IFD T1 centcurs verify well with tasic expectations Lased on 
: physical reasoning and theory. Surface reflecticr and 
7" kottcr refraction patterns are observed as anticipated. Far 
field ncxraal modé prepagaticn can te verified in the ficts 
using sisple dump distance caiculations. In short, all veri- 
ficaticzr nethods attempted, faii to uncover any incorsis- 
tency in IFI performance in a shallow water environment. 
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III. [ABORATOCHY MEASUREBENIS 


a. EACKCECOND 


Ihe najcr attempt at appraising IFD performance invwclved 
comparizg uodel results with lakoratory measurements. The 
Shailcw water envirorsent mcdeled in tae tank is very ideal- 
ized; a ten degree sloping sand bottom with an iscvelccity 
water ccluagn. Althcugh this scenario appears exteresly 
Simplistic, it is ore€ that can ke reasonably modelec ir the 
laboratcry and still approximate conditiors in an actual 
Shallcw water ocean ervironnent. The sethods used tc nodei 
and sé€asure the accustic field are relatively untested. 
Indeed, this attempt at lakoratory modeling was ;erfcrued 
not cnly tc verify IFC predictions, but also to see if the 
€nvircrment could be successfully modeled in the latecratcry. 


E. EXEEBISENTAL DESICN 


A fiberglass coated wocden tank was used. Ibe tank 
is 304 centimeters ir length, 117 centimeters wide ard 95 
centizeters deep. Sand filling the bottom of the tank was 
shaped tc form tke ten deégree slcping bottou, and 
measurements were taken over a range of two meters. Maxizun 
water depth in the tank was 35 centimeters. A 100 kiz 
Source sas placed at mid-channel depth (17.5 centizeters). 
The laycut cf the tark is depicted in Figure 3.1. 


| « 
o 
r=) 
s 
v 
wn 
: a rb 
. => c = 
O O 8 
E 
: a ° 4 
> ° nu s 
- a a 
; - 2 
w“ E HH ; 
j = Oo v risen 
Zz ue + se ete 
; . - es 
: ) ee 
= hed vee 
- fea) a 
My = < ae 
s a) ae 
- ro} nu ; 
. ” Ye) ee 


2 
ae 
- 
. 


| 
te 


* 43 


Fi et a et a te 
ple “° eee i ee 


ete 
hel et Se 


A slope of ten degrees was selected for several 
reascns. Ic begin with, even though a ten degree slcse is 
greater than most ocean pottcn slopes, it is stiil saali 
enough tc ke considered realistic. Pernaps nost imfortact, 
the ten degree slcje was selected because given the 
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frequency limitations (discussed iater in this chapter) and 
range lixitations, this wedge anyle allowed the source tc be i 
placed many (41.6) dump distances rror the apex. A large 


number ci dump distarces was necessary to simulate a distart 
source. 
Ihe bottcm material used in the experiment was #30 ae 


fine crade sand. The grain size ranged fron 0.15 millimeters 
to 0.70 silimeters. The sand was treated wita a techniyue 
used ty Eaek (1984) tc remove air from tne sediment. This 
technique used a high speed jet to agitate the sandyxater 


wey 


mixture to remove tke bubbles and then allowed the sard to 
settle fcr several days before the experiment was initiated. 
i. . Fresh(tar) water served as the medium in the tack. 
Io regove air bukbles, the water was allowed to settle ina 
settling tank before being transferred to the experimental 
tank fer use. The water in the tank was periodically treated 
with chicrire bleach to prevent the growth of biclcgical 
gsaterial. 


- 


a Z. Signal Generating/Receiving Equipment 


The acoustic signal used for the measurements was 
Froduced ky a functicn generatcr, sent through an aunplifier 
and thea transmitted into the water by a directional trars- 
2 ducer rescnznt at 100 kHz. The dimensions of the active face 


* cf the transducer were 7.0 cm in width and 2.0 cm in height. * 
2. These ditensions resulted in an approximate beamwidth (anyie ie 
from the acoustic axis to the first theoretical null) of 
y 11.9 degrees in the horizontal and 46.3 degrees ir the 


er. vertical. Ihe narrcw horizcotal beam minimized reflecticns 
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from tke sidewalls cf the tark, while the wide vertical rean 
allcweéd ccmplete ersonificaticn of the channel in the 
vertical dizensicn. 

Ike signal was received ky an LC-5 omni-directioral 
hydrcghcre, sent threugh an agplifier ard filter, and then : 
displayed cn an oscillosccre and voitmeter. A schenatic 
showing the electronic setur is Shown in Figure 3.2. 

A freguency cf 100 kHz was selected tor tue frac- 
ticai reascns. First, to avoid particle scatterirg bt; the 
sediment, the acoustic wavelength gwust be at least three 
times larger than tle grain size of the sedinent (Andérsor 
and LlieLermanrn, 196€&). The iargest grain size in the sand 
was C.C7 cn, So that a wavelenta of 1.45 cm(100 kEz) was 


sufficiently large erough to be immune to this effect. 
Second, the 100 kHz freguency and the properties oi the sand a " 
Frovide a dump distarce that is smali enouga to alicw the : 
source tc ke positiored many dump distances from the apex. 
Shaging the sand bottom into a ten degree wedge with 
a unbifcrm and smootk interface proved to ke aklecng and 
tedicus precess. Ic facilitate this modeling, wooden 
supperts (two-by-fours) were mounted azony the length of 
koth sides cf the tark. These supports were elevated at cne 
end cf the tank to achieve the required ten degree slope. A 


scraping device was ccnstructed with wooden supports anda 
metal scraping tlade that extended across the width cf the 
tank. This scraping device consisted of wooden suyyorts 
along the top that reached across the tank and cculd be 
fulled alcng the elevated wocden supports on both side of 
the tank. This scrapirg device was pulled along the suprorts 


repeatedly until a szooth slore of ten degrees was sculf- 
tured frca the sand. 

Ecles had tc be drilled into the metal scraping 
Klade tecause when a solid tlade was used in the shallow 
Forticn cf the slope, water trapped behind the tlade was 
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Figure 3.2 Electrenic Eguipment Schematic. 
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f forced ceneath the scraper gouging the smooth bottor. This 
Fottcr rcdeliny techrigue was slow because, once a fas¢ sas 
made cver the bettcz, the water became turbid and it sas ‘ 
then in,cssible to see the rotton. When the bottom was rot 
visifle, it was imygcssible tc see where further saccthaag 
WaS Necessary until the water settled several hours later. 


In addition, as this smoothing process continued a silty 
residue kecame separated frem the sand and settled cut on ‘ 
top cf the sand. Tkis residue would te easily resuspended bor 
and eventually had tc ke removed using a water syphon. 


Cw MEASUREMENT EFEROCELIURES 


Measurement cf tke pressure field within the water sas 
done ty lcwering the receiver in depth at specific ranges of 
interest. The receiving hydrophone was attached to a pair of 
gicrcreters at rigat angles to one another, that was in turn 
rolted tc a board which spanned the width of the tark. Once ae 
the tcard was placed close to a range of interest, cne i, ee 
wicrcseter was used to give fine adjustments in range and — 
the ctker in depth. : 

The measurements were subject to both an accuracy anda : 
Frecisicn error. On a given day, with the water level fixed oe 
and the cress-tank support set at a particular place in the no 
tank, it was possible to fosition the receiving hydrcficne es 
with an accuracy in depth and range of plus or minus 0.06 wey 
centineters (one turr of the micrometer). To prevent the pores 
sand ‘inland' of the apex frem drying out when not taking 
measurements, enough water was added to the tank after ¢€ach 
data run tc keep the sand ccrpletely submerged. The next 
time se€asurements were taken, water had to to te removed oo 
from the tank tc reestablish the beach. Because cf these 


small changes in the water level the horizontal fositicn cf “ 
the Kkeack was subject to a precison error estimated to be Mele 
withir plus or minus cne centineter. eee 
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Ike fitral Gecisicn asscciated with the aeasurerzerts 
centered cr whether tc use a triggered pulse cr a ccrtinucus 
wave (C¥) signal. With a triggered pulse it was possirle to 
distirguish the received signal from interference caused ry 
reflecticnse off the side of the tank. On the other Larc, by 
using a triggered fulse there was a possizility that the 
Fulse length was not long encuga for the acoustic erergy 
associated with paths reflecting off the top and Lcttcr of 
the water cclumn to cvérlap the direct path from the scurce 
to receiver. A CW Signal weuld avoid pctential pulse lenyth 
problems, tut it would be impessible to distinguish Letweer 
the actual signal ang interference. The use of a source wath 
a harrew hcrizontal ream reduced the effects of reflecticrs 
from the side walls, but there was still the possibility of 
interference from side lobes reflected from the sides. 

It was necessary to determine the best zeans tc take 
measurements. This sas done ty takiny measurements at the 
same locaticn with different fulse lenyths to determire the 
teguired rulse length to give consistent results and then 
comparing these results to those obtained with a C& signal. 
The tkird dump distarce (14.4 cm from the apex) and jist 
fast the tenth dump Gistance (50.0 cm from the apex) were 
chosé€n, and measurements were taken for pulsed signals cf 64 
and 2£6 cycles, and acC# signal. These measureterts can be 
seen in Figure 3.3 through Figure 3.5. In the fiytre 
Gepictirg results at the third dump distance, the CW 
measurements are shown as the solid curve, the 64 triyger 
cycle results are disslayed as circular points, and the 256 
trigger cycle results are dericted as triangular pcints. for 
Clarity, the measurements taken just past the tenth dump 
distarce are shown ir two figures. Figure 3.4 compares Ci 
results (sclid curve) with the 64 triyger cycle sigral 
(circular joints). Figuré 3.5 compares CW aeasurenerts 
(solié curve) with the 256 cycle results (circular foints). 
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Figure 3.3 Eulse Length Analysis at 3.0X. 
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Figure 3.5 Eulse Length Analysis at 10.4x. age 
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The pressure amplittde values were normalized Ly dividing 
each [ressure value ty the maximum fressure value ir the 
field. Ike depth values were also normalized by dGiviéing 
each depth ky the maxinun depth. 

The results in Figure 3.3 show gocd agreement Fetweer 
the three curves. In particular, the results icr the Cw 
signal ard the 256 cycle signal are very sinilar. Ike 64 
cycle signal shows tle same behavior as the CW sigral, rut 
differs in magnitude Lelow rid-depth. 

The results at £C.0 cm fren the apex appear much nore 
complicated than at the third dump distance. Figure 3.4 
compares results ackeived with the Cw sigcal ard the 4 
cycle sicnal. From the figure it can be seen tnat the two 
sets cf results shcw good agreement. The measurements 
depicted in Figure 3.£, which ccmpares the 256 cycie sigral 
and the CW signal show foorer agreement, althcugh tke 
general shares cf both curves remain similar. 

Ikese results ixdicate that interference from tne side 
walls is sufficiently small sc that it is possible tc make 
measurezents with a Ck signal to at least ten dump distances 
from the apex. Long pulses could also te used, but the 
difficulty cf makinc voltage geasurements with an cscillo- 
scope ccmpared to reading a voltmeter, dictates tkat 
geasurementse should te made with a CW sagnal. 
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A. IBTECLCUCTION 


IFD acdel jfredictions and laboratory measurements were 
cbtaired fcr compariscn in Eroth a generai and detailed anal- 
ysis. The general analysis compared results every five 
centizeters from the beach to a range or 50.0 ca fires the 
reach. The measurements were spaced to give results at 
a pprceximately each cf the first ten dump distances. The 
detailed analysis ccmrgared results every centimeter frco 2.0 
to 11.0 cm from the ayex. These measurements were taken to 
cbserve the sensitivity of the pressure field tc smal 
changés ir range. 


Be TEE CENERAL ANALYSIS 


Ike general anlysis compared IFD results and lakczratcry 
measurements at appicximately eacn of the first ten dump 
distarces (every five centimeters from the beach). these 
measulezents were oktained by fixing the receiving hydro- 
Fhone with respect tc the crcss board and then moving the 
roard out in increments of five centimeters to measure the 
desired field as a functicn of depth at each range. the 
received CW signal was read on the voltmeter. 

Tke ccafarisens are shown in Figure 4.1 through Figure 
4.10. In the figures, the IFD predictions are displayed as a 
solid line while the experimental measurements are shcwn as 
circles ccnnected by a dashed line. Each pressure value was 
nermalized to unity ty dividing it by the maximum fresstre 
value fcr the respective curve. The depth values were aiso 
normalized ty dividing each depth by the maximum depth in 
the water ccluana. 
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Figure 4.2 Comparison of ReSults at 2.1X. 
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Figure 4.3 Comparison of Results at 3.1%. 


56 


NORMALIZED AMPLITUDE (P/PMAX ) 
0.2 0.4 0.6 0.8 1.0 


Figure 4.4 Comparison of Results at 4.2X. 
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Figure 4.6 Comparison of Results at 6.2%. 
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Figure 4.7 Ccmparison of ReSults at 7.3X. 
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Figure 4.10 Comparison of Results at 10. 4X. 
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: In general, the ;ressure fatterns predictei ry the 1FD 
‘. and those taeasured, koth tbecome more compiicated as the 


| i Tange frem the Leach increases. At all ranges, taere 1S 

z gualitative agreement in the scale of the predicted features 

: : and tke scale of the measured features. Quantitative agree- 

gent is lacking. This agreement in scale but not detail, 

i suggests that the phases amcng the normal anodes predicted by 

. the scdel do not accurately reflect the experinenatal a 
~ situaticr. : 
| Cw. TEE LETAILED ANALYSIS a 


Ike detailed analysis ccmpared IFD values wita exceri- 
mental measurements every centimeter from 3.9 cm to 11.9 ca 
from the teach (0.7X tc 2.3X). ‘Zhe laboratory neasurerents 


were taken by fixing the toard at one Ilccation ard taea bates 
using the wicrometer to adjust the receiver to the desired “ “ 
Tange. Tke received Ck signal was read from the voltmeter. — 
Ikese results are depicted in Figures 4.11 tarough 4.19. ae 
a In tke figures, the IFD predictions are shown as a solid mn 
a ; line and tke experisental measurements are displayed as nee 
= circles. Each jressure value was normalized to unity by a 
-. dividing it by the maximum pressure for the respective a 
7 curve. The depth values were normalized by dividing eaca 
depth ty the maximum Gerth at the particular range. 
5 Tkere is gualitative agreement in the scale of tke kLasic 
ie features fcr all rances, but guantitative agreement is rot 
4 cbserved. Tkte IFD patterns change very little tarcughcut the 
ey analysis, while the measured values chanye nore rcazicly 
- (especially past 2.0%). The results at Z.1X and 2.3x ne 
2 (Figure 4.1€ and Figure 4.19) Show that the pressure field ee 
2 can change fairly significantly over a range as short as cne eee 
~ centigzeter. At those ranges where there is poor agreement a 
- tetween results, there is an indication of  fhase ae 
: 
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Figure 4.12 Comparison of Results at 0. 8X. 
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Figure 4.13 Comparison of Results at 1. 0X. 
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Figure 4.14 Comparisen of Results at 1. 3X. 
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Figure 4.15 Conpariscn of Results at 1. 5X. — 
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interference (Figure 4.18). For distances less than SX, 
theory predicts that for a source at mid-depth, the lcwest 
propagating mode shculd experience interference frem only 
the €vanescent tails of higher modes. Therefore, ir the 
regicr cf the detailed analysis, oue expects smail 
interference effects. The experimental resuits in this 
regicn however, shcw rather Significant interference. the 
interference suggests unsuspected propagating modes are 
preset, er that tke evanescent tails are larger than 
expected. From the trend in the curves it aprears as if the 
Fhase interference is not a factor from 1.5X inward tcward 
the keack (Figures 4.11 througb 4.15). 
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V. COSCLUSIONS/BECOMMENDATIONS 


a. CCNCIUSIONS 


1. Eerformance cf the IFD Model 


This analysis of the IFD acoustic moiel did rot. 
uncover any major failures cf model performance in a Sigpli- 
fied shallcw water environment. But this is not to say that 
tne gscdel consistently and accurately performs in such an 
ccean scenario. Althcuga comparison of IFD predicticns witha 
two ctkher acdels and with simple physicai reasoning did not 
uncover any inconsistencies in performance, the agreement 
rketween IFD values and labcratory measurements is insuffi- 
cient te give complete confidence in the performance cf 
either tke sodel or tke experiszent. 

There is reasonable agreement between the scale of 
the features precicted by the IFD model and the experizert. 
Eoth the model and the mneéasurements show increasing 
complexity as range is increased from the beach. But despite 
the sizilarities, there is poor quantitative agreemert in 
results. Cre possikle cause for the differences may ke that 
the phases of the ncrmal mcdes predicted by the ecdel are 
extreszéely sensitive tc minor irregularities in the share of 
the interface and tke acoustic properties of the [cttcn. 
Frelizinary work (LeSesne, 1984), suggests that the fhase 
relaticnshigrs betweer modes is strongly dependent cn the 
distarce of the source from the apex even at great ranges. 
Consecuently, it agyears that the collective influence cf 
the ncerual sodes is dependent upon careful geometric ccatrol 
cf the experiment. 

Ike detailed analysis around 2X, reveals that large 
Fhase interferences cccur where only one propagating mcde is 


ets . 

. -_"e woe “ hen er er Pr, a en a ae “Le ' 
im i to OS, eat at of ole Fe at at J e . . 
iy ° 


expected. Tnis indicates the existenc2: of anh extéenied 


évanescert tail withir cutoff ci the nigher wodes. Kavéezera 
and Icanreu (1978), noted that the apparent phases cf tae 
evanescert tails are eéxtrenely sensitive to the jietails of 
the ervirtcraent. The results indicate that the tail dces rot 
decay quickly, and its influence is pervasive (Figures 4.2 
and 4.3). From the trend in the curves, it agrears this 
phase interference is not a factor any closer tc the leach 
than 1.5X (Figures 4.11 through 4.15). 


a“ 


Zz. ModelingsMeasurement Erccedures 


The labcratcry measurements represent an initial 
attempt at rodeling ar idealized shallow water 2uvircrment. 
The experimental techri,ues are not without problems. All 
sets cf measurements were repeated and snowed gcod agree- 
ment. Tke reproducitility cf the measurements suggest that 
randca errers have teen minimized. However, Systematic 
€rrors remain which also contribute to the discrefancies 
ketween predicted and measured values. Some of these errcrs 
derive frcm ecuigpment, such as the fact that the size cf the 
hydrcyhone is of the same order of magnitude as the scale cf 
variaticrs in the pressure field. In addition, the envircn- 
ment in the tank may net be sufficiently close to the ideal 
envircngeéent assuzed ry the rodel. 


E. RECCHEBERDATIONS 


Fertrer study and verification of the IFD computer model 
is reccurended. A cc&farison of IFD predictions with ctkter 
models nct analyzed in this study, may uncover the cause of 
the inccnsistency wken compared to measured values. A 
detailed study of the phase interaction of the normal nci¢s, 
althcugh extremely ccmplex, may also offer insight intc the 
IFD ferfcrnance. 
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It appears that the mecdeliny of an idealizeji shallew 
water envircnment in the lakoratery is of vaiye in acoustic 
podel verificaticn. [Lespite the experimental difficulties, 
there was qualitative agreerent in the tasic features 
retween mrcdel predicticns ang laboratory mseasureazents. £Eut 
the €xpéerizent represents only an initial jrobe into tne 
nodelingyméasuregent technigues of a sShailow water envirca- 
nent. Further refinement of these technigues and the exzeri- 
qental eéguipment may result in better guantitacive 
PUP eR ES betwee:. acdel crediciions aia e2steLaiec:. 


acadotieczents. 
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APEENEIZ C 
FOMNING IHE CCNTOUR FLCT ON THE NPS COMPUTER 
A. INIBCLUCTION 


This Agrencix descrites a frocedure for running the [fL 
contcur slet on the NES computer. Detailed instructions for 
Tunnirg the IFD srogram can be found in Jaeger(1983). 


E. CCEYING THE FILE FOR USE 


Cnoce tke TfL values have been generated by the IFD 
Frogranr, aii that is needed to frodauce a II contour flict is 
the fiie FICTS FCKETRAK. This file is shcwn in Appendix E and 
can te copied frog a computer account maintained ry the 
Undersater Acoustics Curriculum. To Jlink with this account 
and cltain a ccry of the file, the user shculd frecceed as 
fcllcsss 


(1) Leg cn terminal. 

(Zz) Enters CP LINK 0160P 181 195 KB .« 

(3) Wken promsted for the read passwerd enters UX . 

(4) Ertexrs acc 19&£ Cc. 

(5) Enter: COFY PICIS FOBIBAN C==A. 
At this fecint the PICIS FOBTRAN file should reside cn the 
user‘s A disk. 
C. kEUBNING THE EROGEAA 


Eefcre runnirg tké progragr the user must obtain the TL 
waluee fxicm the IFD frogram. These TL values pust Le sent to 


415 


Pr i Te i Rae S 


wep 


2S Se Se Sw ee ee em ae ee Se te San Se Me Le Bae ee ee ee an a a 


: a data disk in crder tor the EICOIS program to te arle tec tse a 
: thes. This can te dcre ty flacing a WRITE statemexrt in the mas 
l : IFD freyras that sends the values to a data disk. The sodi- we 

: fied IFD dericted in Appendix A uses this technique ané can ok 
: ke used as a guide. on 
: Cuce tke data disk has teen created, the user smtst ae 
I assign tespcrary disk space (IDISK) to give the pregran 

“ sufficiert room to generate the IL contours. For mere infcr- 


Baticn cn how to assign temporary disk space, see IPS 
Technical Ncte TN-VM-C1 which is availible in the computer 
i consultart's office. om 
With these initial steps ccapleted, all that remains is ice 
to ceapile and run tke program. The prograp can be ccmzpiled 


OFM eee 


rys 

Enter: FCKTGI PLOTS . eae 
. Tae frograg must be run at the TEK618 graphics térgital oe 
under DISSFIA. This can be dene ty: peat 
a Enter: DISSPLA . a 
“ beat 
a The user will then te prempted for the compiled Fertran te 
1 Frogras rame and the iile definitions for the data éisk, nen 
iy tefore tke program will run. 
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APEENRIZ 2 
SOUSCE LEPIH SENSITIVITY ANALYSIS 


ji .i- + the laboratcry experiment, two suppl3 mental sets 
cf measurements were taken to obtain an indicaticr cf the 
Sensitivity of the fressure azplitude to changes in the 
source depth. The first set of measurements was oktained by 
varying the source défth sith the receiver fixed at 1.0X 
(4.8 cm frcem the beach), and lowered in depth to the Ecttca. 
For tke seccnd set of measurerents, The receiver was fixed 
at the third dusp distance (14.4 cm from the beach), and 
geasulements were taken with the source fixed at 5, 7, 9, 
11, 13, 15, 17, 19, 21, 23, 25, 27, 29, and and 31 cs fron 
the strface. 

Ihe results of the first analysis are shown in Figtre 
Dei. In the figure, fressure was normalized by dividirg by 
the gaxigsum pressure and depth was normalized by dividirg by 
the degth cf the water colugn. Although theory predicts only 
cne propagating mode at this distance, the results show 
modal rhase interfererce. 

The results cf the seccnd analysis are shown in Figtre 
0.2 tkrcugh Figure D.€. In the figures, each curve repfre- 
sents a set of measurements taken with the source fixed at a 
specific derth. For convenience, more than one curve is 
shown cn a civen figtre. All fressure values were norgralized 
ky dividing by the tke saxiauzt pressure in the field. the 
depth values were ncraalized ty dividing by the saxiaua 
water degth. 

It yenéral, at :.0X there appears to be sodal rhase 
interference which is influenced Ly source derth. Giver the 
rougbness cf this analysis, the details of this interference 
are ckacure. However, the aodal interference does nct afgrear 
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Figure D.1 Source Sensitivity Analysis at 1.0%. 
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Figure D.2 Neasuresents with Source Depth of 5, 7, and S Ca. 
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Figure D.3 Neasuresents with Source Depth of 11, 13, and 15 Ia. 
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Figure D.4 MSeasurexzents with Source Depth of 17, 19, and 21 Ca. 
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